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Abstract 

The Schroedinger equation for the charmonium and bottomonium states at finite temperature is 
solved by employing an effective temperature dependent potential given by a linear combination of 
the color singlet free and internal energies obtained on the lattice from the Polyakov loop correlation 
functions. The melting temperatures and other properties of the quarkonium states are evaluated. 
The consistency of the potential model approach with the available lattice data on the quarkonium 
temporal correlators and spectral functions is explored. 
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I. INTRODUCTION 



The anomalous suppression of the J /if) production in heavy ion collisions — experimen- 
tally observed in the depletion of the dilepton multiplicity in the region of invariant 
mass corresponding to the J/if> meson — has been proposed long time ago as a possibly 
unambiguous signal of the onset of deconfmement [3f. Indeed, in Ref. js!] it was argued that 
charmonium states can only be produced in the first instants after the nucleus-nucleus colli- 
sion, before the formation of a thermalized Quark-Gluon Plasma (QGP). Then, in their path 
through the deconfined medium, the original cc bound states tend to melt, since the binding 
(colour) Coulomb potential is screened by the large number of colour charges. This, in turn, 
would produce an anomalous (with respect to normal nuclear absorption) drop in the J /if} 
yields. In this picture it was implicitly assumed that, once the charmonium dissociates, the 
heavy quarks hadronize by combining with light quarks only (recombination leading to a 
secondary J /if) production being neglected). 

Since in hadronic collisions a sizable fraction of the measured J/^'s comes from the decay 
of excited charmonium states (bottomonium feed-down contribution being negligible at the 
Super Proton Synchrotron (SPS) conditions) and the latter are expected to dissociate at 
lower temperatures, a mechanism of sequential suppression has been proposed, with the aim 
of reproducing the J /if) suppression pattern as a function of the energy density reached in 
the heavy ion collision (the highest temperatures and energy densities being reached in the 

nnnnn 

most central collisions) [4], |5|, |6|, |7|, ISj . 

In all the above mentioned works the reason for charmonium melting in a hot environment 
was essentially found in the Debye screening of the potential among the two color charges. 
In this regard one immediately faces a few theoretical problems. 

First of all, it is not even clear, at present, whether a temperature-dependent potential 
can describe medium modifications of quarkonia, although there are efforts aiming to de- 
rive a potential at finite temperatures from the underlying QCD (see, e. g., Refs. [9|, Il0|). 
Pragmatically, one can follow a phenomenological approach, building an effective model and 
testing it against quarkonium properties directly calculated in QCD. However, even assum- 
ing the validity of the potential model at finite temperatures, one has to tackle the problem of 
finding the appropriate effective screened potential to insert into the Schroedinger equation, 
as it will be discussed below. 



2 



Indeed, when the mechanism of anomalous J/ip suppression was originally proposed, 
lattice calculations could not provide any quantitative estimate of the Debye screening length 
and the authors had to employ a schematic model for the cc interaction. Nowadays very 
precise lattice calculations of Polyakov line correlators are available, for different numbers 



of light dynamical fermions 
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12. 



Q, Q, Q. 



It has been known for a long time [3] 
that from these correlators it is possible to extract — in hot QCD, as a function of the 
temperature and of the QQ separation — the change in free energy once a QQ pair is placed 
in a thermal bath of gluons and light quarks. 

The free energy obtained in these calculations, taken in the color singlet channel, has been 
used as a temperature-depende nt p otential and inserted into the Schroedinger equation 



in a number of works 
temperatures for al 



13, 



18(. Such a choice turned out to give very low meltint 



the charmonium states: the dissociation temperatures Td = 1.10T C |6| 
and Td = 0.99T C [17| were found for the J/ip, all the other charmonium states melting well 
below T c . 

On the other hand, since the free energy contains an entropy contribution, it was soon 
realized that employing the change in internal (instead of free) energy as an effective QQ 



potential could appear better justified 



19 
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22] . This choice results in a more attrac- 



tive potential, leading to a melting temperature for the J/ip around 1.5 -j- 2T C , the other 
charmonium states [ip' an d Xc) dissociating a bit above T r . These fin ding s appear in agree- 
ment with (quenched) lattice spectral function studies [23, 



24 



25 



27] . at least for what 



concerns the melting temperatures of charmonia. Lattice results are starting to be available 



also for the bottomonium spectral function 



27 



29] and for the unquenched case 



They appear also able to explain the most recent analysis of the Na50 (Pb-Pb), Na60 (In- In) 



3JJ and RHIC (Au-Au) 22j data on J/i[) production, which seems to favour a scenario in 



which only the excited states of charmonium are anomalously suppressed. 



Other choices for the effective potential can be found in the literature 32|, |33|, [34j. In 



Refs. 



32 



in particular, the author argues that the internal energy obtained from lattice 
data through the standard thermodynamical relation contains also a gluon and light quark 
contribution. A subtraction procedure is constructed, leading to an effective QQ potential 
given by a linear combination of the lattice free and internal energies (normalized to the 
case with no heavy color sources), with temperature dependent coefficients obtained from 
the QCD equation of state. 



3 



Results for the quarkonium dissociation temperatures, in potential models based upon 
lattice data, appear to be consistent with independent lattice results from spectral functions 
calculations. Indeed, because of the huge numerical effort required by lattice spectral studies, 
a very fine temperature scan is not available yet. Furthermore, since — as mentioned above 
- a derivation of the potential model at finite temperature from first principles is still 
lacking, it would be desirable to test it against other observables. 

Actually, full information on the fate of the QQ states in thermal equilibrium in a hot en- 
vironment is encoded, depending upon the channel one is considering, in the meson spectral 
functions. The latter can be extracted from the euclidean temporal correlators of mesonic 
currents (measured on the lattice) by inverting an integral transform. Such a task is usually 



achieved with a technique referred to as Maximum Entro py Method (M 



23 



24 



25 



26 



J, 



35IJ36 

1Q. 



361 ] and 



lattice results using this method can be found in Refs. 

Since the reliability of MEM to extract the spectral functions is not fully established 
yet, the check of consistency of the potential model has been mainly devoted to a direct 



comparison with the euclidean correlators in a number of papers |37l . 1381 . 1391 . 140L 14 ll . 142 ] . 
In these works the authors, starting from different screened potentials, calculate in a given 
channel the corresponding charmonium spectral function, which gets contributions both 
from bound states (as long as they are supported by the potential) and (starting from a 
threshold energy) from scattering states. Convoluting the spectral function with a thermal 
kernel, they eventually obtain the charmonium correlator along the imaginary temporal 
direction. Such a quantity can then be compared with the ones measured on the lattice. 



In Refs. 



37 



38 



39 



401 ] the authors point out that this procedure leads to a disagreement 



between lattice and potential model results. Hence they conclude that a study of the QQ 
states in the QGP in terms of screened potentials may not be able to catch the right physics. 



On the other hand, the authors of Ref. 



411 ] propose of curing the discrepancy (at least in the 



pseudoscalar channel) by keeping in the continuum part of the potential model spectrum 
only the resonant contributions, i. e. by subtracting the free gas states. We believe that 
this procedure might be incorrect, since the evaluated correlator has to be compared with 
the lattice ones, which do have a free gas (infinite temperature) limit. 

The comparison between lattice data and potential model calculations is usually done 
by considering a quantity constructed ad hoc, to display the temperature dependence of 
quarkonium properties, namely the ratio between euclidean correlators above and below 
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the critical temperature (see next Section). In a potential model this ratio turns out to be 
very sensitive to the treatment of the continuum, e. g. to the threshold energy. This is 
clearly a problem for the calculation of correlators below the critical temperature, since, for 
instance, at T = the potential model is based on confining potentials and the continuum 
spectrum has to be added by hand. Secondly, lattice calculations appear to be dominated 
by finite size artifacts right in the continuum region, leading to unphysical peaks in the 



30. 



431 ] . Moreover, the asymptotic 



spectral functions, even in the infinite temperature limit 
high energy behavior, <j{uo) ~ u 2 , of the continuum spectral functions is not reproduced by 
lattice calculations, since the finite lattice spacing provides an ultraviolet cutoff. 



In this paper, following Refs. 



32 



331 ]. we extract from the lattice data for Polyakov 



line correlators an effective temperature-dependent QQ potential and we use it in order to 
understand to which extent a comparison with mesonic temporal correlators obtained on 
the lattice can be pursued. In particular, we study the effect of different models for the 
continuum and we try to keep under control the uncertainties introduced by the need of 
calculating correlators below T c by employing the effective potential derived from lattice 
data for T < T c as well. 

We wish to stress how important (both from a theoretical point of view and for its 
experimental consequences) would it be to show that an approach in terms of screened 
potentials is indeed able to catch the main features of the behavior of quarkonia in a hot 
environment. 

Indeed, at present, from finite temperature lattice meson spectral functions one can 
extract important information on the fate of the ground state in the different quantum 
number channels, concerning its position, strength and melting temperature. Unfortunately, 
no reliable information can be extracted so far on the excited states, due to the difficulty 
of disentangling their peaks (if any) from lattice artifacts. Although this is not a severe 
limitation for charmonium studies, where one expects just a few states to (possibly) survive 
after the phase transition (J/ip, Xc and ip', the latter being the only excited state), this is 
no longer true for the bottomonium. 

From the experimental point of view this turns out to be a strong limitation, especially 
in view of the ALICE experiment, where one expects a sizable production of bb pairs. In 
particular the vector channel states will be experimentally accessible through their elec- 
tromagnetic decays into dilepton pairs Q|. Moreover, charmonium yields will also get a 
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feed-down contribution from bb states. The knowledge of the melting temperature of the dif- 
ferent states is then important in order to develop a reliable model of sequential suppression 
to be eventually compared with the experimental data. 

In this regard, it is clearly of importance, in view of quantitative studies of bottomonia 
above T c , to demonstrate that potential model calculations are consistent with the results 
obtained from the lattice spectral functions. 

Although, as we shall see, definite conclusions cannot be drawn yet, our work constitutes 
a promising consistency check of the potential model, given the present status of lattice 
calculations at finite temperatures. Of course, the final aim should be, on the one hand, 
to establish a clean theoretical justification of the potential model and, on the other hand, 
to employ it in the interpretation of the next generation of heavy ion experiments. As a 
general remark, we remind the reader that potential model approaches, like the one adopted 
in this paper, give a static picture of the behavior of charmonium in the QGP, the effects 
of the interaction with the deconfined quarks and gluons being encoded in the screened 
QQ interaction. Alternative pictures are available in the literature, which attempt to pre- 



dict final charmonium yields in nucleus-nuc 



dissociation-recombination processes 
of hadronization 



45 
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eus collisions as arising from the interplay of 



48 



49 



471 ] or in the framework of statistical models 



50J. 



The paper is organized as follows: in Section [TT] we first briefly review the formalism of 
finite temperature quarkonium spectral functions and imaginary time correlators and then 
the procedure followed to extract the effective QQ potential from lattice data; in Section II III 
we compare the outcome of the potential model to results from lattice calculations, both 
for the spectral functions and the euclidean correlators; finally, in Section IIVI we summarize 
and discuss our results. 



II. FORMALISM 

A. Lattice correlators and spectral functions 

The standard object introduced in lattice studies of quarkonium properties at finite tem- 
perature T is the Euclidean time correlator, defined as the thermal expectation value of 
a hadronic current-current correlation function in Euclidean time r for a given mesonic 
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channel H (see, e. g., Ref. [5l|, Chap. 7): 

G H (T,T) = (j H (r)jU0)), (1) 

where jn = Q^hQ and T H = 1, 75, 7^, 7^75. The four vertex operators T H correspond, 
respectively, to the scalar, pseudoscalar, vector and axial-vector mesonic channels, which 
in turn, at zero temperature, correspond to the Xco (xt>o), Vc (Vb), J/^ (T) and Xd (x&i) 
quarkonium states for the cc (66) system, respectively. Although the correlation functions 
can in general be defined for any spatial momentum p, most lattice studies are restricted to 
the case p = and we shall consider, in the following, only this case. 

The correlators of Eq. ([I]) are directly evaluated in lattice QCD, whereas physical observ- 
ables are related to the spectral function an, i- e. to the imaginary part of the real time 
retarded correlators. The two quantities are connected by an integral transform 

POO 

G h (t,T)= duja H (uj,T)K(T,u,T), (2) 
J 

which is regulated by the temperature kernel 

cosher -1/2T)] 
K(T '"' T)= sinh[^/2T] ' (3) 

the energy integration extending over the whole spectrum. 

Extracting the continuous spectral function from the discrete and finite — and usually 

rather limited — set of lattice data is an ill-posed problem. It has been tackled through 



the application of the MEM 



35 



361 ] . which appears to yield promising results, although its 
reliability has yet to be confirmed. 

The MEM analysis has been applied by different groups to the cc system, both below 
and above the critical temperature and mainly in the quenched approximation (see, however, 
Ref. [30] for a study in presence of dynamical fermions). The spectral functions obtained in 
these studies present, schematically, similar features: a well defined peak in correspondence 
to the ground state meson in a channel of given angular momentum (S- or P-wave) up to 
the dissociation temperature; inability to resolve radial excitations; presence in the highest 
part of the spectrum of peaks that are associated to lattice artifacts. Actually, the most 
meaningful information extracted from the lattice spectral functions concerns the existence 
and mass value of the lowest quarkonium states at a given temperature. 

However, the survival, e. g., of the J/^f up to temperatures around 1.5-j-2T c has important 
consequences for the interpretation of heavy ion experiments and it has raised the need for 



a confirmation not suffering from the uncertainties of the MEM analysis. A procedure has 
been devised that makes directly use of the euclidean correlators, by comparing the behavior 
of the correlators above and below T c 23j. Specifically, one introduces the ratio between 
the correlation function Gh{t, T>) of Eq. (J2J) at a temperature T> > T c and the so-called 

reconstructed correlator 

GfffaT^TJ = J dwoH^TjKfru)^), (4) 

calculated using the kernel at T> > T c and the MEM spectral function at some reference 
temperature T< < T c . This procedure should eliminate the trivial temperature dependence 
due to the kernel and differences from one in the ratio should then be ascribed to the tem- 
perature dependence of the spectral function. Furthermore, the MEM result for ^(a;, T<) 
- because of the larger number of lattice sites at low temperature — is more robust. On 
the other hand, at high temperature the correlator is directly measured on the lattice, thus 
avoiding the extraction of the spectral function with the MEM procedure 

Indeed, lattice studies show that the ratio Gh/G t £ c stays around one up to the same 



dissociation temperatures extracted from the MEM ana 



ysis and then it departs from one 



23]. 



following different patterns in the different channels [23 . 

The other approach followed in the literature in the determination of the quarkonia disso- 
ciation temperatures — which is based on the use of effective potentials for the quarkonium 
systems — has reached similar conclusion about the dissociation temperatures, at least when 



making use o 
temperature 



effective potentials extracted from lattice data on QQ free energies at finite 
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33|. 



In a potential model approach one solves the Schroedinger equation for a given potential, 
thus getting not only the spectrum but also the wave functions of the QQ pair. It is then 
fairly straightforward to get the corresponding spectral function as the imaginary part of 
the QQ propagator, namely 

a H (u,T) = -ImG H (u) 

7T 

= ^2\{0\jH\n)f 5(u - K) 

n 

= F U(" - M n ) + 6(u - s )Fl U} _ SQ , (5) 

n 

where n represents all the relevant quantum numbers (in the second line the sum over n 
actually stands for both the sum over discrete states and the integration over the continuum) 
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and E n = so + e n , so being the continuum threshold and e n the eigenvalue associated to the 
state |n). In the last line we have separated the discrete and continuum contributions, 
introducing the couplings Fn, n = (0|j_ff|n) and F H e = (0\j H \e) and the bound state masses 
M n (a normalization of the discrete and continuum states as (n\n') = 5 nn r and (e|e') = 
5(e — e'), respectively, is understood)^ 



Following Ref. 52] (see also Ref. 



39J) one can express the couplings in terms of the wave 



function at the origin for the ^-states (pseudoscalar and vector channels), 

F 2 PS = ^\R(0)\ 2 and F 2 = ^\R(0)\ 2 , (6) 

and in terms of the first derivative of the wave function at the origin for the P-states (scalar 
and axial- vector channels), 

F 2 s = ^- 2 \R'(0)\* and F\ = ^|i?'(0)| 2 , (7) 

N c being the number of colors and m the quark mass. 

By inserting the spectral function of Eq. ([5]) into Eq. ([2]), one gets the potential model 
expression for the Euclidean correlators: 

POO 

G H (T,T) = y2FH,nK(r,M n ,T)+ de F 2 . e K(r, e + s ,T). (8) 

Some care is necessary in dealing with the continuum part of the spectral function. Indeed, 
in Eq. ([8]) the integration over the excitation energy can reach pretty high values, especially 
for small r's, probing regions where the non-relativistic dispersion relation is no longer valid. 
Actually, the asymptotic relativistic spectral function uh(uj) that one gets in perturbative 
QCD is known to be proportional to uj 2 in every channel, whereas it is easy to check that the 
non-relativistic solution of the Schroedinger equation in the free case goes like, e. g., u 1 ^ 2 for 
the 5-states. The difference is due to the different phase space and it can be accounted for 
by simply renormalizing the wave function with a phase space factor 53(. In the calculations 



discussed below we will be using such wave functions, renormalized to account for relativistic 
kinematical effects. 

Note that in potential model calculations of the spectral functions, a mixed representation 
has been sometimes employed, with a continuum part approximated using t he p erturbative 



QCD expression with a threshold extracted from the effective potential 



37 



39j: 



a H (u,T) = y^FlJiu - M n ) + —u 2 9(u - s )f H {w, s ), (9) 



where the function /# is given by leading order perturbative calculations 54|, 

f H (uJ,so)= an + bn 1 ^ \U- 



55 



56j as 



(10) 



At leading order, the coefficients (a#, &#) are (1, —1), (1, 0), (2, 1) and (2, —3), for the scalar, 
pseudoscalar, vector and axial-vector channels, respectively. This form for the spectral 
function is affected by an obvious inconsistency between the bound and the continuum 
parts of the spectrum. Nevertheless, in the following we make use also of this expression in 
order to test the model dependence of our results. 



B. Effective potentials from lattice data 



In Ref. 20] we have provided a unified parametrization of the temperature and separation 



dependence o 



' the lattice data for the color-singlet QQ free energy F\ in the case of quenched 



ll|, 2-flavor [l4j and 3-flavor [l5( QCD. 

In the following we employ the parametrization obtained there for the case Nj = - 
since most lattice calculations of Euclidean correlators are in quenched approximation - 
but we shall also consider Nt = 2, since in this case lattice data for the free energy are 

n 

available to us also at temperatures 0.76 < T/T c < 1 [14| and one can use them to study 
the behavior of the correlation functions through the phase transition. We refer the reader 
to Ref. [20[ for a detailed description of the fitting procedure. 

Although in past work the free energy Fi has been directly used as an input for the 
QQ potential energy, more recently it has been recognized that the QQ internal energy 
Ui provides a more appropriate candidate. Since the two quantities are connected by the 
well-known relation 

F = U-TS, (11) 

once a suitable parametrization of the temperature dependence for the free energy has been 
obtained one can subtract the entropy contribution and get the color-singlet internal energy 
as 

Vl = -T>™ , 12 ) 



The latter has been employed in Ref. [20] as the effective QQ potential in the Schroedinger 
equation, getting good agreement with the charmonium dissociation temperatures extracted 
from lattice calculations of the spectral functions. 
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T/T 



FIG. 1: 

N f = 2 



n 

'he ratio a(T) = 3p/e as a function of T for Nf = [57( (T > T c , dashed line) and 



58|] (T > 0.76T C , solid line). 



However, the use of the internal energy (1121) is not fully satisfactory yet. From the 



theoretical point of view, as pointed out in Ref. [32| and thoroughly discussed in Ref. 33], the 
internal energy of Eq. ffl2|) contains the genuine QQ potential energy, but also a contribution 
due to the variation of the gluon (and light quarks, for Nf ^ 0) internal energy in the presence 
of the QQ pair. 

Also the comparison with independent lattice results presents some shortcomings. In 
fact, although it is quite successful in providing dissociation temperatures in agreement with 
direct lattice estimates, its asymptotic (r — > oo) value — employed to define the threshold 
energy entering into the continuum part of the spectral function — turns out to have a 
sharp temperature dependence around the critical temperature. This in turn is reflected, as 
we shall see in the next Section, into a temperature dependence of the correlation functions 
which is not observed on the lattice. 

In this work we follow the procedure described in Ref. 32|] in order to separate the true QQ 



internal energy, , from the gluon and light quark contribution. The effective potential 
entering into the Schroedinger equation is defined as (F 1 and U\ are the same quantities 
considered in Eq. (fl2l) ): 



Vi(r, T) = U? Q (r, T) - Uf Q (r -> oo, T) 

= f F (T)[F 1 (r,T)-F 1 (r^oo,T)] + f u (T)[U 1 {r,T)-U 1 (r^oo,T)}, (13) 
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where 



f F (T) = — (14a) 
6 + a(i J 

A/CO = -^2=, (14b) 
3 + a(i J 

and 

= (15) 
p and e being the pressure and energy density of a homogeneous system of quarks and gluons, 



respectively. These thermodynamical quantities have been obtained in quenched [57J and 
2-flavor [58] QCD as a function of temperature and we display in Fig. [1] the resulting ratio 
(j!5p in the range of temperatures we are interested in. Note that the weight functions fp 
and fu turn out to vary in the range 3/4 < fp < 1 and < fu < 1/4, respectively, a fact 
which explains why the QQ potential is closer (but not identical) to the free energy F\ than 
to the internal energy U\. 

The asymptotic value of the QQ internal energy, {r — > oo,T), is used to define 
the continuum threshold, Sq(T) = 2m + U® {r — > oo,T) and we compare it in Fig. [2] 
to the asymptotic value of the internal energy as obtained from Eq. (fl2!) . i. e. without 
subtracting the gluon and light quark contribution. The sharp variation of U\ around T c 
- which, as discussed in the next Section, leads to correlation functions in contrast with 
lattice calculations — is evident, whereas displays a much smoother transition. 

The effective potential based upon is less attractive than the one based upon U\ 
(but still more attractive than the one based upon i*\), a fact which gives rise to slightly 
lower dissociation temperatures. Values for the latter in accord with the results from lat- 
tice spectral function studies can be obtained by using quark masses slightly higher than 
the values of the Particle Data Group listing 59(. For instance, in Refs. 33[ dissociation 
temperatures compatible with the lattice phenomenology have been found using a charm 
quark mass of 1.4 4- 1.5 GeV. Note that the lattice free energies employed to parametrize 
the effective QQ potential have been obtained for infinitely heavy quark mass and that the 
lattice calculations of correlation functions use a variety of quark masses, usually chosen in 
order to reproduce the mass of the charmonium (bottomonium) ground state. 

Here we solve the Schroedinger equation using an effective quark mass defined as rh(T) = 
m + (T)/2, which, beside improving the binding of the QQ pair, can be naturally 
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N=0 



0.75 



1.25 1.5 
T/T 



1.75 



FIG. 2: The asymptotic values of the lattice free energy [14] (dotted), the total internal energy of 
Eq. (fT2|) (dashed) and the QQ internal energy (solid) in 2-flavor and quenched QCD. 



interpreted as a thermal mass. Using this definition the continuum threshold is given as 
Sq(T) = 2m(T). Our choice implies that the same thermal mass is employed both in the 
kinetic term of the effective hamiltonian and for the energy of the heavy quark at rest 
in the plasma. In principle, nothing prevents these two masses from assuming different 
values. While the threshold mass can be extracted from the available lattice data, in order 
to evaluate the thermal correction to the kinetic operator one should know the momentum 
dependence of the heavy quark self-energy in the QGP. Lacking the above information, here 
we choose to ignore such a distinction. Note that in a non relativistic QED plasma a moving 
test charge receives a positive thermal correction to its mass, which is slowly decreasing 

i this rigorous result 



39| that the thermal 



with the temperature [60]. Our choice is in qualitative agreement wit 
obtained in another contest. On the other hand, it has been argued 
mass should not be used in the Schroedinger equation, since in the bound state the heavy 
quarks should not feel the effect of the medium. Yet, the quantitative differences between the 
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TABLE I: Spontaneous dissociation temperatures (in units of T c ) of the lowest cc and bb states. 
The values in square brackets have been taken from Ref. [20J (see text); the ones labeled with (*) 
have been obtained with a potential extrapolated beyond the temperature range of the lattice data 
(T < 2T C for Nf = 2). 
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< 1 


1.00 


1.00 
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< 1 


< 1 


0.98 


0.99 




m b = 4.3 GeV 


m b = 4.7 GeV 


m b = 4.3 GeV 


m b = 4.7 GeV 




2.96 [4.5] 


3.18 


3.9(*) [6.7(*)] 


4.4(*) 


Xb 


1.13 [1.55] 


1.15 


1.15 [1.63] 


1.17 




1.12 [1.40] 


1.14 


1.13 [1.43] 


1.15 



two approaches are quite moderate, well below other uncertainties present in the problem. 
III. RESULTS 

In this Section we use the effective potentials of Eq. (1131) to calculate spectral and corre- 
lation functions by solving the Schroedinger equation 

4> £ (r,T) = e(T)i; e (r,T), (16) 

for the QQ eigenvalues e(T) and the eigenf unctions i^eir, T). The dissociation temperatures 
for the lowest cc and bb states are shown in Tab.[H Since the potential is spin-independent the 
two 5-states, pseudoscalar and vector (or P-states, scalar and axial-vector) are degenerate. 

With respect to the previous findings with the full internal energy [20] one now observes a 
reduction of the dissociation temperatures: in the cc channel the ground state melts around 
w 1.5 -T- 1.6 T c , whereas the excited states already around the critical temperature; in the bb 
channel the ground state melts above 3T C and the excited states around ~ 1.1 -v- 1.2 T c . For 
comparison, in the Table we also report (in square brackets) the dissociation temperatures 
obtained in Ref. using the full internal energy, in the cases where the same quark mass 
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V 2 

■— + V 1 (r,T) 
m 



had been employed (bottomonium for m& = 4.3 GeV). 

In Figs. [3] and 0] we show the mass, M = 2m + Ui Q {r — ► oo, T) + e(T), of the lowest 
S-w&ve and P-wave states, respectively, as a function of temperature. Apart from a narrow 
range of temperatures around T c , the variation of the masses with T is quite moderate, the 
difference between the maximum and minimum values being around 150 -j- 200 MeV of the 
iS-states and ~ 100 MeV for the P-states. This, for instance, amounts to roughly 5 4- 6% 
(2%) of the J/ip (T) mass. 

A slight decrease of the r] c mass in the r ang e of temperatures explored here has also been 
observed in some lattice spectral studies |6JJ. Note that a more pronounced softening of 
the J/ty peak (here degenerate with rj c ) as the temperature increases has been attributed 
by some authors to the presence in the vector channel of a transport contribution to the 
spectral function, not resolved by the MEM procedure, but nevertheless tending to move 
strength towards lower energies 27]. 

Just above T c one observes a slight peak in the masses for the case of 2-flavor QCD. Note 
that this is the range of temperatures where the T-dependence of the parameters employed 
to fit the lattice free energies is stronger [^J, so that it might just signal the inadequacy of 
the parametrization at the critical temperature. It is however curious that the same behavior 
does not appear for Nf = 0, where the same parametrization has been employed. 

Note that the smooth temperature dependence of the masses of the QQ states is strictly 
related to the smooth T-dependence of the asymptotic QQ internal energy, (oo,T), 
obtained through the procedure described in Section III BL By employing, in the definition 
of the bound state mass, the total internal energy Z7i(oo, T), one would get much higher (by 
a few GeV) masses around the critical temperature (see Fig. |2J). 

Other observables can be easily calculated in the potential model. As an example, we 
show in Fig. [5], for the cc system, the square of the S'-wave radial wave function, P(0) 2 , 
and the square of the first derivative of the P-wave radial wave function, P'(0) 2 , calculated 
in the origin. This quantities can be related to the quarkonium decay widths for different 
processes (see, e. g., Ref. [51]), such as the leptonic decay rate of a neutral vector meson 



63, 



64[ | or x decay into pairs of pseudoscalar or vector mesons 65]. 



As one can see from the figure, the values of R(0) and, particularly, of R'(0) are quite sta- 
ble below T c and very close to the values at T = 0. Above T c the P-wave states have melted, 
whereas R(0) for the 5-state drops almost linearly up to the dissociation temperature. 
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FIG. 3: Mass as a function of temperature of the lowest 5-wave cc (upper panel) and bb (lower 
panel) states obtained from the solution of the Schroedinger equation. The upper curves are for 
a quark mass m c = 1.6 GeV (m& = 4.7 GeV), the lower ones for a quark mass m c = 1.4 GeV 
{nib = 4.3 GeV). The arrows point to the value of the mass at T = 0. 
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FIG. 4: As in Fig. [3l but for the lowest P-wave states. 
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FIG. 5: Squared value in the origin, for the cc system (m c = 1.6 GeV), of the S-wave radial wave 
function (left panel) and of the first derivative of the P-wave radial wave function (right panel), as 
a function of the temperature, for quenched and unquenched QCD. The arrows point to the values 
at T = 0, calculated with the Cornell potential of Ref. 62]. 



A. Spectral functions 

Here we would like to compare the spectral functions calculated in the potential model to 
the ones extracted from the lattice euclidean correlators using MEM. To make easier a direct 
comparison of the spectral functions it is convenient to modify the potential model expression 
of Eq. (j3j) in order to accommodate a width for the bound states. Indeed, the width of the 
bound states at finite T in the lattice spectral functions contains a contribution due to 
the (possible) thermal broadening of the states and a contribution due to the statistical 
uncertainties, the latter apparently being dominant above T c because of the limited set of 
data points [2^\ - 

A width can be easily included by substituting in Eq. (J5j) the delta function in the bound 
state sector with a Lorentz distribution: 

where T n represents the width of the state n. 

A word of caution is in order about the calculations at T < T c — where we do not 
deal with a quark-gluon plasma, but with a still confined system. Although the use of a 
potential model for bound states is well assessed, at least at T = 0, the interpretation of 
the continuum spectrum is much more questionable. Here we shall use both the spectral 
function calculated in perturbative QCD, as usually done in the literature, and the one 
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FIG. 6: The J/ijj spectral function divided by uo 2 as a function of lo for Nf = 2 and Nj = at 
several temperatures below and above T c ; in the left panels the continuum part of the spectrum 
comes from the solution of the Schroedinger equation, in the right panels from the perturbative 
expression. 
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coming from the potential model, in order to get a feeling about the model dependence of 
the results. We would like however to stress that the aim of the model is to describe the 
physics of quarkonia at T > T c and that calculations in the confined phase are done in 
order to discuss the lattice results for the correlation functions, which are mainly available 
in the literature as ratios of functions above and below T c (see Sections III Al and IIII Bp . 



In potential model calcu 



at zero temperature 
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ations, the reconstructed correlators have been usually considered 
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421 ] . adding to the bound state contribution — calculated 
with a phenomenological confining potential — the perturbative QCD continuum. This 
choice obviously poses another problem of consistency (besides the one, already mentioned, 
of the consistency between the bound and the continuum spectra), since above T c one is 
using a lattice-derived potential and at T = a phenomenological (fitted to the data) one. 
Incidentally, the phenomenological parameters of the T = potential adopted in the above 
works (i.e., the coefficients of the Coulomb and linear terms) turn out to be quite different 

n 

with respect to the ones employed in normalizing the free energies at short distances |14| . 
The use, also below T c , of lattice-derived potentials should give one a feeling about these 
model dependencies. 

In Fig. [6] we display the cc spectral function for the S-states for quenched and unquenched 
QCD, both below and above T c (since we neglect the hyperfine spin-spin interaction the 
J/ijj and r] c states are degenerate and the corresponding spectral functions differ for trivial 
factors). In the left hand panels we employ — for the continuum part of the spectral function 
- the solution of the Schroedinger equation (Eq. ( 1TT1) . see the discussion in Section Hi Al) . 
whereas in the right hand panels we employ the perturbative expression (see Eq. ([9])). In all 
cases a width T n = 100 MeV has been used for the discrete states. Here and in the following 
calculations we employ m c = 1.6 GeV and m;, = 4.7 GeV for the quark masss. 

In this figure one can note a marked peak of constant strength and almost fixed position 
below T c , corresponding to the J/ip state, and a second smaller peak, corresponding to 
its first radial excitation. Above T c only the ground state peak survives and its position 
moderately moves (consistently with the mild mass shift displayed in Fig. [3]), whereas its 
strength is gradually decreasing as one approaches the dissociation temperature. There is 
no significant difference between Nf = and Nf = 2. 

Above T c the bound state peak is clearly visible especially in the right hand panels, 
where the perturbative continuum spectral function has been employed. The continuum 
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FIG. 7: The J/V> spectral function divided by u? as a function of oj for Nf = 2 and T = 1.4 T c 
(dotted line). The bound state contribution (solid line) is also shown separately. 
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FIG. 8: Bound state contribution to the pseudoscalar spectral function divided by J 1 for Nf = 2 
and Nf = at various temperatures. A width r„ = 700 MeV has been used. 
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FIG. 9: The Xc spectral function divided by J 1 as a function of u for Nf = 2 at several temperatures 
below and above T c ; in the left panel the continuum part of the spectrum comes from the solution 
of the Schroedinger equation, in the right panel from the perturbative expression. 
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contribution calculated from the Schroedinger equation, on the other hand, shows a resonant 
part when the bound state energy is approaching zero and this resonant contribution tends 
to dominate over the bound state peak (see Fig. [TJ where, as an example, the spectral 
function at T = 1.4 T c is magnified around the threshold energy, separating the bound state 
and continuum contributions). 

In order to show more clearly the evolution of the bound state above T c — and to make 
easier a comparison with the lattice results — we display in Fig. [8] only the bound state 
part of the S- wave spectral function using a larger effective width, comparable to the one 
of the lattice charmonium spectral functions calculated in Refs. 23, 27] (Nf = 0) and [30} 
(Nf = 2). As one can see the pattern as a function of the temperature for the position of 
the bound state appears to be similar to the lattice results. 

The temperature dependence of the strength of the peak, instead, seems to differ from 
the lattice outcome. In Refs. 23|, |27j, the 5*- wave spectral functions have also been extracted 



at various temperatures using the same number of data points, in order to compare results 
affected by the same statistical uncertainties: at least for the pseudoscalar state (the vector 
channel being more uncertain) it has been found that the spectral function maintains its 
strength up to 1.5 T c . This result is, of course, affected by large errors and also depends upon 
the reliability of the MEM procedure. On the other hand, the potential model calculation 
yields a practically constant strength up to T c and then a practically linear decrease, as it 
could be inferred also from Fig. [5j 

In Fig. [9] we show the spectral function for the P-states in the case of unquenched QCD, 
using, for illustration, a width of 100 MeV (again, we neglect hyperfine splitting of scalar 
and axial- vector states). The continuum contribution in the left panel comes from solving 
the Schroedinger equation, in the right panel from the perturbative calculation. Also in 
this channel one can note, below T c , a marked peak of constant strength and approximately 
fixed position, corresponding to the \c state. No bound states are present above T c : the 
peak appearing in the left panel at T = 1.05T C is actually a resonance in the continuum 
contribution. Also the narrow peaks visible at higher energy for T < T c are embedded in 
the continuum. The situation of the P-wave lattice spectral functions is less clear, since 



this channel is much harder to study on the lattice. All the studies [23l . 1271 . |30j of the cc 
system observe a strong modification of the spectral functions above T c and this is usually 
interpreted as a signature of dissociation of the P-wave states. However, above T c some 
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FIG. 10: Bound state contribution to the spectral function divided by uj 2 for bb states at temper- 
atures below and above T c . Nf = 2 in the left panels and Nf = in the right panels. The widths 



are 650 MeV (T) and 1 GeV (xb) 



28, 



29] 



strength is observed at very low energies, especially in the scalar channel, well below the 
expected continuum threshold and there is no physical explanation for these findings. 

Finally, in Fig. [TD]we display the bound state part of the S- and P-wave spectral functions 
for bb states, in the range of temperatures studied in Refs. [2^] and 29]- The T state appears 
to be stable both in position and strength over the whole range of temperatures, whereas 
the Xb state gets strongly modified approaching the dissociation temperature (see Tabled]). 
The stability in the position of the T and the dissociation of the Xb at a temperature 
T = 1.1-i- 1.2T C is in accord with the results from lattice calculations; the latter, however, 
are still preliminary and one cannot draw from them any conclusion about the evolution of 
the strength with T. We do not show the continuum part of the spectrum, since it is very 
similar to the charmonium case. 

To summarize our findings for the QQ spectral functions, we can state that the bound 
state part of the potential model spectrum appears to be consistent with the lattice results 
as far as the evolution with T of the position of the cc and bb states is concerned. On the 
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other hand, both the evolution of the bound state strength and the shape of the continuum 
contribution — either perturbative or calculated in the potential model — show major 
differences from the lattice outcome. 



B. Correlation functions 

As we have already mentioned above, the reliability of the MEM procedure is still under 
discussion. In order to reveal the temperature dependence of the spectral function, it has 
been proposed to consider the ratio between the euclidean correlator at a given tempera- 
ture and the reconstructed correlator, which contains the spectral function at a reference 
temperature (see Section Hi Aj) . 



In lattice calculations, this has been done using as reference temperature both some finite 
T < T c [23] and T = [27]. In potential model calculations, the spectral function at T = 
- modeled as in Eq. ([9]) by using a phenomenological potential for the bound states and 
the perturbative QCD contribution for the continuum states — has been employed so far as 



a reference 
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421 ] . However, such a spectral function is not only affected by the inconsis- 
tency between the bound and the continuum states, but it also strongly depends upon the 
threshold energy sq, which has to be treated as a parameter. Moreover, employing, for the 
continuum contribution in the confined phase, the perturbative QCD spectral function is at 
least problematic. Here we try to remove the first two sources of uncertainties by employing 
the spectral functions calculated from the lattice-generated potential at T < T c . Of course, 
the problem of interpretation of the continuum spectrum is still there, but, as we shall see 
below, it is probably inessential. 

In Fig. [11] we display the ratio Gh / G r ^ c for the pseudoscalar and scalar charmonium 
states at different temperatures, using the reconstructed correlator at T = 0.75T C . Both 
the continuum spectra generated by solving the Schroedinger equation (left panels) and by 
using the perturbative expressions (right panels) are employed. 

Let us first analyze the ratio at T < T c . Here it turns out to be always very close to 
one, not only for the case shown in the figure, but over the whole range of temperatures 



available to us 
on the lattice 



0.75 < T/T c < 1). It strongly resembles the ratio of the correlators measured 



23(| : even the small decrease at large r's in the scalar channel seems to be 



correctly reproduced. Note that this outcome simply reflects the temperature independence 
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FIG. 11: Ratio Gu/G r ff for the pseudoscalar (upper panels) and scalar (lower panels) charmonium 
states (Nf = 2) at different temperatures, using the reconstructed correlator at T = 0.75T C , as a 
function of the euclidean time r. In the left panels the continuum part of the spectrum comes from 
the solution of the Schroedinger equation, in the right panels from the perturbative expression. 



of the spectral functions below T c , as it was apparent in Fig. [6j This evolution with the 
temperature in our model stems from the mild dependence upon T of both the effective 
potential and the threshold energy, which, in turn, is due to the procedure of subtraction of 
the gluon and light quark contributions from the QQ internal energy (see Section fll Bp . The 
same quantities calculated using directly the internal QQ energy as the effective potential 
would yield a strong T dependence already below T c (see, e. g., Fig. |2J). 

Above T c , the ratio in the pseudoscalar channel looks qualitatively similar to the one 



measured on the lattice 



23 



27] 



at variance with the result of Refs. 
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42|] 



but 



there are important differences. In the lattice results, the ratio remains close to one up to 
the dissociation temperature and then it gradually decreases; here, the departure from one 
occurs as soon as the temperature grows beyond T c and reflects the rapid change in the 
threshold energy visible in Fig. [2j After that the smooth change in the threshold energy and 
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the reduction of the bound state strength as T increases give rise to the moderate decrease 
of Gh I G T £ C observed in Fig. [TTJ 

In the scalar channel the situation is rather different: the ratio calculated using the 
perturbative continuum grows above one and then drops down, whereas in the case of 
the potential model continuum it is always lower then one. Both results are at variance 
with the lattice ones, where the ratio is uniformly growing with r. The reason for the 
different behavior employing either the perturbative or the "interacting" continuum can be 
reconducted to the presence, in the latter case, of a strong resonance just above the threshold 
at T < T c (see Fig.EJ. 

As one can see, differences in the outcome from the different models employed here and in 



other works 



39 



421 ] can in general be understood in terms of differences in the treatment of 



the continuum contribution to the spectral functions (its form and/or the threshold energy). 
Hence, the ratio Gh/G^ c appears to be a sensitive observable in discriminating among 
different models. However, so far no model calculation has been able to reproduce the 
lattice results for every channel. 

In order to understand the reason for this discrepancy, we display in Fig. [12] the spectral 
functions (divided by a; 2 ) in the pseudoscalar and scalar channels at a fixed temperature. 
In the same plots we also show (dot-dashed lines) the factors u 2 K(t,uj,T) by which the 
spectral function a H (u,T)/u 2 should be multiplied in order to get the integrand yielding 
the euclidean correlator of Eq. ([2j). The factors u 2 K(t,u,T) are shown for two values of r, 
namely a small one — where the ratio Gh/G^ c goes to one both in lattice measurements 
and in model calculations — and a large one — where the effects of temperature on the 
ratio are stronger. 

As one can see from the figure, at small r's the weighting factor u 2 K is a rather smooth 
function covering a wide range of energies and the corresponding correlator gives, roughly 
speaking, a measure of the total strength associated to the spectral function. On the other 
hand, for large values of r the weighting factor uj 2 K is a rapidly dropping function of the 
energy, selecting only the low energy part of the spectral function (this is, of course, a nice 
feature for discriminating among different models). 

The discrepancy with the lattice results can be understood by comparing, e. g., Fig. [12 
with Fig. 17 (pseudoscalar case) and Fig. 19 (scalar case) of Ref. [27J. In the pseudoscalar 
case, at T = 2T C no bound state is present in our potential model calculation. The lattice 
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FIG. 12: Spectral functions (divided by lo 2 ) for Nf = in the pseudoscalar (T = 2T C ) and scalar 
(T = 1.16T C ) channels as a function of the energy; both the potential model (solid line) and the 
perturbative QCD (dashed line) results are displayed. The dot-dashed lines represent the factor 
uj 2 K(t, u, T) (in arbitrary units) for two values of r. 



estimates at this temperature show a strong dependence upon the default model used to 
extract the spectral function, yielding either a significant change with respect to T = (see 
Fig. 17 of Ref. [2?]]) or almost no temperature dependence. In any case, the small reduction 
(~ 5%) observed in the ratio of correlators is due to (moderate) differences between the 
T = and T = 2T C spectral functions at e nerg ies < 10 GeV, which are dominated by the 
lowest energy peak visible in Fig. 17 of Ref. 27J. Our potential model calculation (solid line 
in Fig. [T2l) also presents a peak (a resonance in the continuum spectrum) — at variance with 
the perturbative QCD spectrum (dashed line) — and this explains the similar behavior of 
Gh/G^ in our calculation and on the lattice. 

In the scalar case, on the other hand, there are no resonances in the potential model 
calculation and no strength is present at low energy (< 4 GeV) — whereas some strength 
appears at low energy in the lattice result (Fig. 19 of Ref. 271 ]) — and this explains the 
different behavior of Gh/G v ^ c for this channel in our calculation and on the lattice. 

We would like to stress that this explanation for the discrepancy between potential model 
and lattice correlators does not rely upon the MEM-based lattice spectral functions, which 
have been used only to illustrate our argument. For instance, the exact lattice spectral 
function in the scalar channel at T = 1.16 T c may not have the form of the MEM-based 
ones (Fig. 11 of Ref. 27J]), but it must have, in the low energy region, more strength than 
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the T = spectral function in order to yield a correlator twice as large as the reconstructed 
one. 

This analysis — together with the direct comparison of the lattice and potential model 
spectral functions done in the previous Section — shows that the discrepancies between the 
two approaches are mostly located in the continuum part of the QQ spectrum. Indeed, in 
this energy domain, the lattice spectral functions (which provide useful information on the 
properties of the ground state) suffer strong limitations, due to the inability of resolving 
excited/resonant states, the presence of unphysical peaks related to lattice artifacts and the 
bad asymptotic high energy behavior. Note, in particular, that the presence of unphysical 
peaks is not related to the MEM procedure, since they appear also in the infinite temperature 



limit 



43]. Hence the ratio Gh/G t ^ c may not be a good candidate to check the consistency 



of the two approaches (potential models vs. lattice studies) . 



IV. SUMMARY AND CONCLUSIONS 



We have given a comprehensive treatment of quarkonia at finite temperature, within the 
framework of a potential model. We have constructed an effective QQ potential as a linear 
combination of the finite temperature QQ internal and free energies — in order to separate 
the genuine QQ energy from the gluon and light quark contributions — and we have used 
as input the QQ free energies obtained on the lattice from Polyakov loop correlators. 

The effective potential that we have obtained yields dissociation temperatures for the cc 
and bb systems that are in agreement with estimates based on independent lattice studies 
of euclidean correlators and spectral functions. Specifically, we have found that in the 
charmonium system only the ground S'-wave states {J/ip and r] c ) survive up to ~ 1.5-^1.6 T c , 
all the other states melting around the critical temperature or below; in the bottomonium 
system, again only the ground S'-wave states (T and rjb) survive up to temperatures > 3 T c , 
while Xb and T' melt around w 1.1 t 1.2 T c . 

In order to gain deeper insight into the comparison of results based upon the potential 
model and lattice studies, we have also calculated the cc and bb spectral functions over 
a wide range of temperatures. To test the model dependence of the results, two models 
for the continuum spectrum have been employed: one using the interacting solution of the 
Schroedinger equation and one using the perturbative QCD spectrum. The evolution with 
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the temperature of the mass of the bound states has been found in good agreement with the 
lattice results, while the T-dependence of the strength of the bound states agrees only below 
T c . Above T c the lattice-based spectral functions show little temperature dependence of the 
bound state strength, whereas in the potential model — where the strength is proportional 
to the coupling Fjj — the latter drops almost linearly. On the other hand, the continuum 
spectrum in the two approaches shows no resemblance at any temperature. 

We have also calculated in the potential model the QQ euclidean correlators, since these 
quantities are directly measured on the lattice and are not affected by the uncertainties 
inherent to the procedure of extraction of the lattice spectral functions. We have found good 
agreement with the lattice results in the pseudoscalar channel and not in the scalar one: we 
have shown that the agreement or disagreement between potential model and lattice in the 
correlators is not driven by the bound state contributions, but by the continuum spectrum. 
Since the latter is known to be strongly affected by artifacts due to the finite size of the 
lattice, the ratio of euclidean correlators, Gh/G t ^ c , does not appear to be appropriate for a 
test of potential models vs. lattice calculations. 

At the present stage, the only substantial conflict between potential model and lattice 
predictions for quarkonium properties at finite temperature seems to be in the T-dependence 
of the bound state strength from the critical temperature up to the melting point, i. e. a 
progressive drop in the potential model and — presumably, since lattice results are available 
only for a few temperatures — a rapid change at the dissociation point in the spectral studies 
on the lattice. Lattice results for the QQ spectral functions above T c , however, are still 
affected by uncertainties related to the MEM procedure and by large statistical errors and 
new measurements with better statistics will probably put this outcome on firmer grounds. 
On the other side, one has to gain deeper insight into the connection of the potential model 
with the underlying QCD at finite temperature, not only to provide a solid and rigorous basis 
for the model, but also in order to extend its applicability (e. g, by including an imaginary 
part to describe the thermal broadening of the quarkonium states). 
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